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Despite the fact that modified theories of gravity, in particular the f(R) gravity models have 
attracted much attention in the last years, the problem of the energy localization in the framework 
of these models has not been addressed. In the present work the concept of energy-momentum 
complexes is presented in this context. We generalize the Landau-Lifshitz prescription of calculating 
the energy-momentum complex to the framework of f(R) gravity. As an important special case, 
we explicitly calculate the energy-momentum complex for the Schwarzschild-de Sitter metric for a 
general f(R) theory as well as for a number of specific, popular choices of f(R)- 
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tin ■ It has been almost a century since the birth of General Relativity and there are still problems that remain unsolved. 

The energy-momentum localization is one of them which till today is treated as a vexed problem. Much attention has 
CN| ' been devoted for this problematic issue. Einstein was the first who tried to solve it by introducing the methodology 
of energy momentum pseudotensors. He presented the first such prescription [l| and after that a plethora of different 
energy-momentum prescriptions were proposed 0, [H, 0, IE S 0] • All these prescriptions were restricted to compute 
the energy as well as the momenta distributions in quasi-Cartesian coordinates. M0ller was the first to present an 
energy-momentum prescription which could be utilized in any coordinate system [§[ . 

The idea of energy-momentum pseudotensors was gravely criticized for several reasons 1,113 ED, El (actually one of 
the drawbacks was the aforesaid use of quasi-Cartesian coordinates which was solved by Miller's prescription) . Firstly, 
' although a symmetric and locally conserved object, its nature is nontensorial and thus its physical interpretation 



I. INTRODUCTION 



seemed obscure 13(. Secondly, different energy- momentum complexes could yield different energy distributions for 
\^ " the same gravitational background [bH . flU ]. Thirdly, energy-momentum complexes were local objects while there 
f*"*. , was commonly believed that the proper energy-momentum of the gravitational field was only total, i.e. it cannot be 
■ localized [l6j . For a long period of time the idea of energy-momentum pseudotensors was relinquished. 

The approach of energy- momentum pseudotensors for the thor ny pro blem of en ergy -momentum localization was 
(Nj rejuvenated in 1990 by Virbhadra and collaborators [13, El, [H, US El 1H, S3, 0, SI, (2(f . Since then, numerous works 
t-H ' have been performed on computing the energy and momenta distributions of different gravitational backgrounds using 
several energy-momentum prescriptions (for a recent list of references see 27]). In 1996 Aguirregabiria, Chamorro 
and Virbhadra [II] showed that five different energy-momentum complexes yield the same energy distribution for any 
J> , Kerr-Schild class metric. Additionally, their results were identical with the results of Penrose 29] and Tod [1(| using 
? ■ the notion of quasi-local mass. Many attempts since then have been performed to give new definitions of quasilocal 
energy in General Relativity [12, H3 HH HH, Hll- Considerable efforts have also been performed in constructing 
superenergy tensors [13] • Motivated by the works of Bel [11, [H, H(| and independently of Robinson [4l|, many 



investigations have been carried out in this field [42j, |43|, |4J, |45|, |4f 

In 1999 Chang, Nester and Chen [47j proved that every energy-momentum complex is associated with a Hamiltonian 
boundary term. Therefore, the energy-momentum complexes can be considered as quasi-local, boundary condition 
dependent conserved quantities. Finally, it should be pointed out that though a long way has been trodden the 
solution to the problem of energy-momentum localization in the framework of General Relativity is way ahead. 
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Another challenge to the development of physical theory of gravitation is that, the plethora of observational data 
collected recently indicates that our universe is undergoing an accelerated expansion. Motivated by this observational 
evidence, we have been in a long hunt for the explanation for this speed-up. Till today, three possible reasons have 
been presented. Two of them, namely the cosmological constant and the quintessence field, are developed in the 
framework of General Relativity. The third one is developed in the framework of alternative theories of gravity. 
In particular, the simplest among the aforesaid models are that in which Einstcin-Hilbcrt action is modified by an 
additional term. 

Modified theories of gravity, especially the f(R) gravity models that replace the Einstein-Hilbert action of General 
Relativity (henceforth abbreviated to GR) with an arbitrary function of the curvature scalar, have been extensively 
studied in recent years (see e.g. [H, |H, [5(3, [5l|, [H, HH, HE ISll an d references therein). The challenges in 
constructing viable models in the light of cosmological constraints (see e.g. [EE [El, [5§| and references therein), 
instabilities [IE [Ml I&3 , solar system constraints (see e.g. [H, H3, [6f| [66j] and references therein) and evolution large 
scale perturbations [67l 1681 [69| are now known. The Solar System constraints are a major obstacle to most theories 
[zE [ZllllS, 111 but they can be completely removed by certain types of models [58|, [66j, |74j, [75j, [76j . 

It is widely known that when a new theory is introduced, it is expected this new theory to successfully answer 
all already-solved (in the framework of the old theory) problems. Moreover, it is anticipated that this new theory 
will be able to address, alleviate, and finally solve problems that the existing old theory cannot. Following this line 
of thought, we address here, to our knowledge, for the very first time the problematic issue of energy-momentum 
localization in the context of f(R) gravity models which as was mentioned above intend to replace GR. We take first 
steps in this direction and consider energy- momentum complexes within f(R) gravity models. 

The remainder of the paper is as follows. In Section II, we present the basic equations and formalism of f(R) gravity. 
In Sections III and IV, the Landau-Lifshitz energy-momentum complex and the Schwarzschild-de Sitter (henceforth 
abbreviated as SdS) metric, or the SdS black hole background, are reviewed. In Section V, we extend the concept of 
the Landau-Lifshitz energy- momentum complex into the framework of f(R) theories. As a special case, we compute 
the energy-momentum complex of the SdS metric for a number of commonly considered f(R) theories. In the final 
section we summarize the results and present our conclusions. 



The action for f(R) gravity is 



II. f{R) GRAVITY FORMALISM 



S= I (fx^l-j— f(R) + C m ), (I) 



16ttG" 



where the standard Einstein-Hilbert action is replaced by a general function of scalar curvature f{R). The corre- 
sponding field equations (in the metric approach) are found by varying with respect to the metric and read 
as 

F{R)R I1V - \f{R)g^ - V»V V F{R) + 9ltv DF(R) = 8^GT™ (2) 

where T™ is the standard minimally coupled stress-energy tensor and F(R) = df /dR. In contrast to the standard 
Einstein's equations from the Einstein-Hilbert action, the field equations are now of higher order in derivatives. 
Contracting the field equations gives 

F(R)R- 2f(R) +2BF(R) = 8nG(p - 3p) (3) 

where we have assumed that we can describe the stress-energy tensor with a perfect fluid. From the contracted 
equation it is clear that in vacuum, any constant scalar curvature metric with R = Rq is a solution of the contracted 
equation as long as F(Rq)Rq — 2J(Rq). In general, the whole set of field equations is solved exactly by the SdS 
metric l3| (for a more recent work on spherically symmetric solutions of modified field equations in f(R) gravity see 
also (Zl) 

ds 2 = Bdt 2 - B^dr 2 - r 2 d6 2 - r 2 sin 2 ^ 2 (4) 

with 

2M , r 2 

B(r =1 A-. 5 

r 3 

The scalar curvature for this metric is Rq = — 4A. Hence any f(R) theory, including the standard General Relativity, 
satisfying the constant curvature condition F(Rq)Rq = 2/(i?o) has the SdS (black hole) metric as an exact solution. 
We will return to this important special case in a later section. 



III. THE LANDAU-LIFSHITZ ENERGY-MOMENTUM COMPLEX 



In general, the energy-momentum complex t^ v (henceforth abbreviated as EMC) carries coordinate dependent 
information on the energy content of the gravitational and matter fields. It sums up the energies of the matter fields 
through the energy-momentum tensor T^ v (henceforth abbreviated as EM), and that of the gravitational field through 
the energy-momentum pseudotensor V 1V (henceforth abbreviated as EMPT), which depends on the coordinate system 
used to describe the system. The EMPT cannot be defined uniquely and a number of suggestions, with different 
mathematical properties, exist. All of them lead to conserved quantities of the gravitational theory. 

The most straightforward conserved quantity is the integrated EMC over the three-dimensional space integral 

Eemc = I d 3 xT 00 (6) 

JB(0,r) 

which represents both the energy of the gravitational field and that of matter inside the coordinate volume -B(0, r). 
In the case of a black hole, it consists of two parts: the black hole mass M and the energy stored in the gravitational 
field t 0Q , therefore 

E E mc=M+ [ d A x^g~t m . (7) 

JB(0,r) 

In the construction of Landau and Lifshitz [2], one looks for an object, rf" va , that is antisymmetric in its indices 
since then d v d a rf va = due to the covariant continuity equation, which in a locally Minkowskian coordinate system, 
simplifies to 

d v T» v = 0. (8) 

Hence 

= d a rf va . (9) 

Einstein's equations give 

K 2 

where k 2 = 8ttG (we have set c = I), and in the locally Minkowskian coordinates Ricci tensor and Ricci scalar reads 
as 



= —{R^ - -Rg^) (10) 



-dadpg^s - d 1 d s g a p) (11) 
and R = d a dPg a p-g al3 ng a p 

= (^V -9 aP 9nWv9c.fi- (12) 

Using these expressions in Eq. (fTOl) . one can rewrite the EM tensor as 0] 

T" v = d a j]^ a (13) 

where 

^ = ^^^[{-9)^-^)]- (14) 
In a locally Minkowskian coordinate system dag^ — and hence one can define 

(-g)T» v ee d a h^ a ee d a dpH» votf} (15) 
where two new tensors, so-called superpotentials, 

h^ a = {-g)^ va = ^71-^ [(-9) isTti* - 9 m 9 vP )] (16) 
= ^[{-g){g^g aP -g^g 1 " 3 )] (17) 



have been defined. 

In a general coordinate system, Eq. [|15p is no longer valid and one defines a new object t^ v such that 

[ _ g){T »» +tn = __. (18) 

The new object, namely the EMPT , is straightforwardly computed in a general coordinate system employing Eq. 
(fT8|) since T^ v can be expressed in terms of the geometric quantities by using the Einstein's equation, i.e. Eq. (jlOp . 
and h^ va is given in Eq. fTB]) . 

Carrying out this somewhat lengthy but routine exercise, one obtains Q 

_ _l (fore pw _ pe pu _ pe pw \( [ia vf3 _ fiv af3\ 
l LL~ 2^|l Zi a/3 i ew 1 au 1 (3e iffi^uJlS i> 9 9) 

_ l _ n fia f3a/yi' ytu r pi/ pw _ pi/ ru pi/ ru ^ 

+rt ea '(r^r^-r^r^)}. (19) 

The Landau-Lifshitz EMC, i.e. r^, can now be evaluated either as a sum of the EM and the EMPT, namely 



r£l = (-g){tZ" L + Tn (20) 



or directly using Eq. (|18[) which is now written as 

nu dflT 



•ll 



dx a 



(21) 



IV. THE SCHWARZSCHILD-DE SITTER METRIC 

In GR, the empty space solution outside a static spherically symmetric mass distribution in a universe with a 
cosmological constant is the SdS metric, or the SdS black hole metric. In spherically symmetric coordinates, it reads 
outside the mass distribution as (in units where G = 1) given in eqs. ([H and ([5]). M is the total mass and A is 
the cosmological constant. Due to the fact that some EMPTs are calculated in cartesian coordinates, we need to 
re-express the SdS black hole metric in cartesian terms. The metric (jlj then reads as 

= B d t^ B -y + / + 2 z2 dx i- x2+ 2 B y + 2 z2 d y> 

x z + y z + z z x z + y z + z z 

x 2 +y 2 + B- X z 2 , 2 (S" 1 - l)2yz , 

- t. ^ ^ — dz - - Trdydz 

x A + y z + z z x A + y z + z z 

(B- 1 - l)2xy , , (B- 1 - l)2xz , , 

- ^-5 ^-z-clxcly-^ -L_^dx d z. 22 

ar + y z + z z x z + y z + z z 

For this metric, but not for the metric in spherically symmetric coordinates, the determinant reads as g = — 1. This 
feature is used regularly in the following as outer factors like \f—g equal unity. 

The EMPTs corresponding to the SdS solution are straightforwardly calculable (in a cartesian coordinate system). 
Working in a space where r > and thus EMPTs and EMCs coincide, we obtain for the 00-component of the 
Landau-Lifshitz EMPT 

00 2 36M 2 + 12MAr 3 + Ar 4 (Ar 2 - 9) 

\9)t 



k 2 r 2 (6M + r(Ar 2 - 3)) 2 
2 -9Ar 4 + (6M + Ar 3 ) 2 
k 2 r 2 (Ar 3 - 3r + 6M) 2 



(23) 



For comparison, one can easily repeat this exercise for the other well-know EMPTs: we find that the Weinberg 
EMPT is equal to that of Landau and Lifshitz, the Einstein and Tolman EMPT is 

tl = -4 (24) 



and the M0ller EMPT is 



2A 

*8 = --y. (25) 



It is evident that the EMPT of the SdS metric strongly depends on the chosen construction, even though EMPTs 
derived in different prescriptions sometimes can be identical. Furthermore, their behavior far from the mass source 
is wildly different. Therefore, by looking only the functional form of a single EMC, it is obscure what is the precise 
physical interpretation of it, even though it represents a conserved quantity. However, we emphasize, that the different 
EMCs differ by a divergence term related to the boundary conditions of the physical situation [47| . 

V. EMPT IN f(R) THEORIES OF GRAVITY 

Among the different EMPTs studied in the literature, Landau-Lifshitz's and Weinberg's prescriptions appear to 
be most straightforwardly suitable for extending into f(R) gravity theories. Here we consider extending the Landau- 
Lifshitz's prescription and leave others for future work. 

In f(R) theories, one can write the field equations as 



T»u = ^ ^f(R)R„u - \g^f(R) - D tl D v f(R) + g IMU nf(R)\ . 



(26) 



Like in GR, the covariant continuity equation holds, i.e. D^T^ V = [82|], suggesting that one should write the RHS 
of Eq. (|26p as a divergence of an object antisymmetric in its indices, i.e. in a form d a b} lva . 

Following Landau's and Lifshitz's prescription and considering a locally Minkowskian coordinate system at a given 
point we obtain 

k 2 T^ = f'(R)R^ - \g^f(R) + (g^g^ - g» a g^)d a dpf(R) 

= f{R)G» v + ±<r(f'(R)R - f(R)) + a [(g^g^ - g^g^)d f{R)} 

= d a [f{R)^rT a + (g^g a0 ~ g» a g v0 )f"{R)d R] - d a f{R)^rT a + ±g^(f'(R)R - f(R)) 

= d a [f{R)Jrr« + (g^g afl - g^g^)f"{R)d p R] + ±g» v (f'(R)R - f(R)) (27) 

where rj^ a is that defined in ([Ft]) . It is noteworthy that the term d a f (R)n 2r q^ va vanishes in locally Minkowskian 
coordinate system, because rf va is linear in the first derivatives of the metrics. We are partially able to write the RHS 
of the field equations, namely Eq. (|27[) . as a divergence. The remaining term, absent in GR, remains problematic 
in a general case without a clear method which would enable us to write it as a four divergence. We can, however, 
proceed in an important special case where the scalar curvature is a constant, R = i?o- The SdS black hole metric 
belongs to such a class of metrics. 

A. The Landau-Lifshitz -energy momentum complex for a metric with constant scalar curvature 

For a metric with a constant scalar curvature, R = Rq, Eq. (127|) simplifies to 

1 

1 

6k 2 
f 

6k 2 



TV = a Q [/'(i? )^ l/Q ] + 77^^(/'( i? o)i?o-/(i?o)) 

d a [f{W va + -hW v x a - 9 m x v ){f '(Ro)Ro - f(Ro))} (28) 



so that the generalized Landau-Lifshitz superpotential takes the form 



hf™ = f{,W a + ^(g^x a - g^x") [f(R )R - f(R )} ■ (29) 



The EMPT t^ v in a general coordinate system defined in the Landau-Lifshitz prescription can now be read out from 
the expression for the EMC (remember, that g = — 1) 

t"" = T" v + t" v = d a h^ a . (30) 
Hence the generalized Landau-Lifshitz EMC reads as 

= f'(RoKl + [f(Ro)Ro - f(Ro)} d a (g^x a - g^x v ) , (31) 
where t££ is the Landau-Lifshitz EMC evaluated in the framework of GR (see Eq. (T2T~|) ) . The 00-component reads as 

r 00 - f'(Ro)T° L ° L + [f(Ro)Ro - f(Ro)] d a (g 00 x a - g 0a x°) 

= f'(Ro)T° L ° L + -^(f(Ro)Ro - f(Ro)) {d l9 m x* + 3. 9 00 ) . (32) 

Eq. (f3"Tj) is a general formula valid for any f(R) theory when the studied metric has constant scalar curvature. The 
standard GR result is recovered when f(R) = R. 

B. Energy-momentum complex of the SdS metric of some f(R) models 
Using Eq. (j3"Tj) we can compute the EMC of the SdS metric in a general f(R) theory: 
r 00 = r(R )4 L + -^(f(R )R -f(R ))(rB'(r)+3B(r)) 



r^(f'(Ro)Ro 
2 -9Ar 4 + (6A/ + Ar 3 ) 2 



= f(Ro)rl° L + — (/'(/2o)i?o - f(Ro)) 3 — r 2 

ok- 2 \ r 3 



K 2 r 2(Ar 3 - 3r + 6M) 2 



f'(Ro) + ^(f(Ro)Ro f(Ro)) (a - ^ - ^^ 2 ) 



(33) 



This result is valid for any f(R) theory that has the SdS metric as a vacuum solution i.e. any theory which satisfies 
the vacuum equation f'(Ro)Ro — 2/(i?o) = 0. Again note that when f(R) = R we recover the standard, i.e. GR, 
form of Landau-Lifshitz EMPT, Eq. (j2"3"l) . as expected. 

An important special case encompassing popular choices of f(R) is a generic action function 



f(R) = R— (-1)"^ + (-l) m - l bR m , (34) 

where n and m are positive integers and a, b any real numbers. This form of function f(R) is widely used in 
cosmological context. In this case the generalized Landau-Lifshitz EMC takes the form 



0-(l+2n) A-n 

_oo _ z i4 - 



r [(12M - 9r + 5Ar 3 )(a(l + n) + b(m - l)(4A)" l+n )] 

(35) 



9r 2 K 2 

[-9Ar 4 + (6M + Ar 3 ) 2 ] [an + (4A) n (4A + 6m(4A)™)] 
9 A(6M-3r + Ar 3 ) 2 

For the form of /(-R) considered above, i.e. Eq. (|3"4"]) , and recalling that for the SdS metric we have Rq — — 4A, 
constant curvature condition can be written as 

(4A)" +1 = a(n + 2) + b(m - 2)(4A) m+n . (36) 

In the special case where m — 2 or b — we get 

(4A)™ +1 

a = ~ — ~TT- (37) 
n + 2 



Note, that not all type (j3"4"|) models are cosmologically viable. It is known, that those vacuum solutions with i?o such 
that f"(Ro) > (note our sign convention) are inherently unstable [6(| and therefore not suitable for cosmological 
model. 

Particularly often used model of f(R) theory of gravity is 

f(R) =R-!L- t R 2 i (38) 

which has a stable vacuum whenever e > 1/(3\/3m 2 )- The 00-component of the corresponding generalized Landau- 
Lifshitz EMC for this model is written as 



1 f, , As r , t1 -9Ar 4 + 6M + Ar 3 2 /x 4 + i? 2 - 2ei? 3 n 

T ^—\r(eRl-2^) 12M - 9r + 5r 3 A -36^ 777^77 o a ^ 21 f> 39 

2 K 2 i? [ J i? (6M - 3r + Ar 3 ) 2 J v ' 

which for the special case of the SdS black hole metric with the cosmological vacuum Rq = —\^3fj, 2 , reduces to 



r oo = 

18r 2 K 2 i? 



r00 = 2 + 373^ f 5 4 3 _ _ 72 [l92^M 2 + 48, 2 Mr 3 - 3 _6,V + ^V]l 

18V3r 2 K 2 1 4 p v ' (24M- 12r + VV^ 3 ) 2 J 



(40) 



Another cosmologically interesting /(-R) gravity model includes also logarithmic dependence on curvature. Thus it 
reads 

f(R) = R+ {-l^^cRT - din f W , (41) 
where its parameters are related to the cosmological constant by the constant curvature condition written now as 



4A 

(4A) m c + 2A + rfln I — 



d+ (4A) m cm = 2 , , 

1 k 

For this model the corresponding 00-component of the Landau-Lifshitz EMC is of the form 

, (d + cm(4A) m + 4A)(-9Ar 4 + (6M + r 3 A) 2 ) 



(42) 



r 



1 { 



18K 2 r 2 L A (6Af - 3r + Ar 3 )' 

+ r (12M - 9r + 5r 3 A) + c(m - l)(4A) m - din f^JJ }• 



(43) 



In any case the generalized Landau-Lifshitz EMPT of a /(i?) model differs crucially from the GR Landau-Lifshitz 
EMPT for the SdS metric. Taking into account the constant curvature condition we can write 

r m = f'(RoK° L + ^f(Ro) (3 - ^ - ^r 2 ) , (44) 

which coincides with the GR Landau-Lifshitz EMPT only if /(i?o) = and f'(Ro) = 1 (implying, due to the constant 
curvature condition, that there is no cosmological constant). This special case, while possible, is not a general property 
of physically meaningful f(R) models indicating that in a general f(R) model the Landau-Lifshitz EMPT will be 
non-trivially related to the corresponding EMPT in GR. For example, it is clear that at large r the two EMPTs have 
different asymptotic limits with t^l ~ f in GR and tll ~ f 2 in a general f(R) model. 



VI. CONCLUSIONS AND DISCUSSION 



The problem of energy localization has been one of the first problems that was treated after the onset of GR. 
Although a number of scientists endeavored to solve it, the energy localization remains a vexed and unsolved problem 
till to date. In this work, motivated by the recent interest in constructing extended models of gravity and in particular 
f(R) gravity models that replace the standard Einstein-Hilbert action of GR, we have introduced for the very first 
time, to our knowledge, the energy localization problem in the framework of f(R) theories of gravity. In particular, 
we have extended the concept of energy-momentum complex in the prescription of Landau-Lifshitz. Although we are 
unable to formulate a completely general expression for the EMC valid for all theories and metrics, we can proceed 



in an important special case where the scalar curvature of the considered metric is constant. In this case, we have 
presented a general formula for the Landau-Lifshitz energy- momentum complex for a general f(R) theory. We find 
that the general relativity result is generalized to encompass an additional term. 

Metrics satisfying the requirement of constant scalar curvature include the Schwarzschild-de Sitter metric, which 
e.g. describes the space-time around spherically symmetric objects in a universe with a cosmological constant. We 
have computed the generalized Landau-Lifshitz EMC for a general f(R) theory that accepts the SdS metric as a 
solution as wells as for a number of f(R) commonly considered in the literature. We find that the GR result is 
generalized by the presence of additional term. The new term is non-trivial as it has a different dependence on the 
coordinate r than the term arising from the GR part. 

It is more than obvious that further study is needed, e.g. other EMC's and their interpretation, i.e. corresponding 
physical boundary conditions, in f(R) models need to be considered. A particularly interesting and a potentially 
fruitful direction to follow in the future is to consider the problem of energy localization in Weinberg's formulation. 
For a more general EMC covering also the non-constant curvature case, a construction of a new type of EMC may 
be a more direct way to proceed as generalization of the Landau-Lifshitz EMC is challenging. The calculations of the 
integrated constants of motion in different models and systems is another example of a relevant open question. Wc 
hope to address these issues in future work. 
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